Universal scaling behavior of interacting nanomagnets 
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In this Letter we report how the influence of the interparticle dipolar interactions, computed by 
the volume sample concentration c, can be described in a general way if scaled by the dimensionless 
parameter Co = 2K/M%. This scaling co-parameter is derived from the energy equation governing 
the system by using simple arguments, and its meaning is analogous to the normalization of the 
influence of the external magnetic field H by the anisotropy field of the particles Ha = 2K/Ms- 
We use a Monte Carlo technique to test the usefulness of this Co parameter for interpreting the 
controversial evolution of the blocking temperature (Tb) vs. c reported in the literature, finding 
that the Ts(c)-curves of different particle types follow the same trend if scaling c/cq. 



Nanosized magnetic materials constitute the 
working body of several nanotechnological applica- 
tions as the widely-used magnetic data storage [l[ 
and the promising nanomagnetic logic Q , which rely 
on a very precise control of the magnetic entities and 
of their relative interaction. The problem is that for 
these closely-arranged systems the dominant inter- 
action is the magnetic dipole-dipole coupling, which 
lacks an analytical solution and exhibits a complex 
interplay with the other more relevant energies in 
the system (anisotropy, Zeeman). Finding a precise 
and easy description of its features arises, therefore, 
as a crucial task to address in the current nanominia- 
turization race towards ever-denser concentrated de- 
signs |3j-|6|. 

Such a purpose has often been pursued in the ba- 
sis of the Neel model Q ■ Dipolar interactions were 
investigated in relation with the anisotropy energy 
barrier of single-domain superparamagnetic (SPM) 
systems and its effect on the characteristic block- 
ing temperature, Tb [8|. The interparticle coupling 
strength was varied by changing the sample concen- 
tration of the system, c. The reasons to choose this 
method are its ease to carry out the experimental 
measurements and to correlate them with the theo- 
retical background [9(. 

Despite the intense research devoted to this field 
in the recent years, no concluding results have been 
yet achieved, finding in the literature very different 
Tg(c) trends for a priori very similar samples. These 
trends range from, i) a regular increase IICM ; ii) a 
rapid increase followed by a saturation [ll[; iii) a 
non-monotonic behaviour 12j; iv) a decrease 131 ]; 
etc. Although analytical models have been devel- 
oped that describe reasonably well specific cases, a 
general model able to account for all the different 
shapes has not been reported so far (see Ref. [Hj]). 

In an effort to shed light on this complex scenario, 
we present in this work an alternative approach: by 



simply rewriting the energies governing the system 
in an adequate way, the influence of the dipolar in- 
teractions in the magnetic behavior of the system 
(introduced as proportional to c) is described in a 
general way by the unitless parameter Co = 2K/Mg. 
In order to demonstrate the universality of this scal- 
ing parameter, we build up our arguments in the 
context of the usual single-domain particles SPM 
frame. 

The SPM model represents the particles by their 
large magnetic supermoment and also by their mag- 
netic anisotropy K — Khi, with K the anisotropy 
constant and hi the easy axis direction. The mag- 
netic supermoment results from the coherent rota- 
tion of the inner atomic moments and is in first ap- 
proximation proportional to the particle volume Vi 
as fti = MsVii^, with Ms the saturation magneti- 
zation and e w the unitary vector that describes its 
orientation. 

The energy per particle for a non-interacting 
system is given by the anisotropy (Ea) and Zeeman 
(Ez) energies as 
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Both energies determine the magnetic evolution 
of the system, being more important their relative 
weight than their absolute values. This is easily em- 
phasized if rewriting Eq. ([1]) as 
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where ifi and Qi are the angles that forms with 
respect to H and K, respectively. Likewise, the 
anisotropy field is defined as Ha = 'i.K/Mg, which 
accounts for the relative importance between Zee- 
man and anisotropy energies. This is highlighted in 
the field dependence of T B oc (1 - H/H A ) 3/2 [3. 



This expression illustrates that the influence of H 
on the particle system does not depend on the abso- 
lute H- value but rather on the ratio H/Ha- Thus, 
particles with different values of K and Ms may 
still follow a similar Tg (H ) shape depending on their 
characteristic Ha value, i.e. Tb = Tb(H/Ha). 

We follow the same line of reasoning to investigate 
the influence of the dipolar interaction. In this case, 
for an iV-particle system the energy per particle is 
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where fy is the vector connecting particles i and 
j, and the summation extends to the neighboring 
particles. Our procedure is to rewrite Eq. (|3|) also 
in terms of H/Ha, trying to find also an easy de- 
pendence for the dipolar energy term. For such a 
purpose, we assume i) a monodisperse system; ii) 
the particles are located into a cubic box of side L, 
so that we may use normalized units ay = r^Nz /L; 
iii) the dimensionless sample concentration c results 
in c = Eili v i/ l3 = NV/L 3 . Thus, the reduced 
energy e^ reads now 
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where we have introduced Co = 2K/Mg. This Cq- 
parameter weights the importance of the dipolar 
interaction energy as c/cq, analogously as H/Ha 
weights the importance of the Zeeman energy. 
Therefore, we can conclude that same that H/Ha 
provides a general behaviour for the magnetic field- 
dependence of the system, it is possible to define also 
a general behaviour for the influence of the dipolar 
interaction, given by c/cq. This result demonstrates 
the existence of an inner scale of energies in the ther- 
momagnetic properties of those nanostructured sys- 
tems as a function of the Ha and Co parameters. 

The physical meaning of cq can be interpreted as 
an estimation of the relative importance between the 
anisotropy and dipolar energies, i.e. Co = 2K/Mg oc 
Ea/E d . It is worth to note that the 2K/M S ra- 
tio had been noticed in the literature as a relevant 
magnitude concerning interacting SPM systems [IB} • 
However, neither its meaning as a general scaling 
factor nor its key-role to solve controversial results 
in the literature (as the Tb vs. c evolution) had been 
reported up to now. 



The next step is to demonstrate that the above 
assertion on Co as a scaling factor for the influence 
of magnetic dipolar interactions can effectively solve 
controversial results in the literature. With this pur- 
pose in mind, we decided to tackle the Tg vs. c prob- 
lem based on its central role amongst theoretical- 
nanoparticle problems as to dipolar interaction con- 
cerns. 

Our procedure consists on comparing the Tb(c) 
data for different types of particles, in order to study 
the different trends reported in the literature. We 
represent the particle types by their correspond- 
ing co-values, and examine their magnetic behaviour 
both as a function of c and c/cq. Since experimen- 
tal data is usually subjected to uncontrolled sec- 
ondary effects (polidispersity jT^], aggregation [13]) 
that might mask the influence of the dipolar inter- 
actions, we have used a Monte Carlo technique that 
allowed us to perfectly control the characteristics of 
the sample and hence to be sure about the nature of 
the different Tg data. The computational procedure 
is the same as described in Ref . 12| , and the results 



are presented in the usual reduced temperature units 
t = ksT/2KV. Tb is roughly evaluated as the max- 
imum of the zero field cooling (ZFC) curves under a 
low applied field H = 100 Oe, as usual, and hence to 
obtain the Tb vs. c data we simulated ZFC processes 
for a systematic variation of Ac = 2.5%. 

To represent different types of particles we char- 
acterized them by their K and Ms values, and asso- 
ciated each type to the corresponding c . To select 
the K and Ms values (and hence the c cases) we de- 
cided to keep a common value of K and vary Ms- By 
performing this way we can expect that in the non- 
interacting limit case Tb will be the same for the 
different particle types, since Tb ~ KV/2bks (i.e. 
independent of Ms), so that we can study the evolu- 
tion of Tb with c from this common non-interacting 
point. It is worth to note that although this is an 
ideal assumption, it is in fact possible to design dif- 
ferent materials so that they have a common (or very 
similar) K value but different Ms- Also, for compu- 
tational purposes it is easier to keep the same value 
of K (included in the reduced temperature units) 
and vary only Ms- The K value is relatively low so 
that dipolar interactions influence easily the system, 
using K = 1.5 x 10 5 erg/cm 3 . The Ms values consid- 
ered and the corresponding Ha and cq parameters 
are shown in Table U together with H/Ha ratio cor- 
respondent to H — 100 Oe. 

The Ms values were chosen so that they repre- 
sent realistic physical values, and the equivalent cq 
cases cover a wide range of Ea vs. Ed, so that cq 
varies between the ^-dominating case (co = 4), 
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TABLE I. M s {emu/cm 3 ) 
co, HA(Oe), and hwo Oe ■ 
1.5 x 10 s erg/cm 3 . 



data and the corresponding 
= 100 Oe/H A values for K = 



M s 



H A h 



100 Oe 



274 4.00 1095 
316 3.00 949 
775 



387 2.00 

548 1.00 

775 0.50 

949 0.33 

1095 0.25 



548 
387 
316 
274 



0.09 
0.10 
0.13 
0.18 
0.26 
0.32 
0.36 



to the i?£)-dominating one (cq = 1/4), covering as 
well intermediate cases. In Fig. [T] we show some 
ZFC curves for specific field and concentration con- 
ditions in order to illustrate the equivalent role of Co 
as a scaling parameter for the sample concentration 
analogous to the Ha value for the magnetic field. 
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FIG. 1. (Color online) ZFC curves for different field and 
interaction conditions, for some Co values: H — 100 Oe, 
c = 0.0 in (a); H/H A = 0.18, c/c = 0.0 in (b); H/H A = 
0.18, c = 0.10 in (c); H/H A = 0.18, c/c = 0.10 in (d). 

Fig. [IJa) displays some ZFC curves for the non- 
interacting case under the same absolute field H = 
100 Oe. In this case, the curves differ from each 
other and also the maximum of the curves; however, 
if using the same reduced field H/Ha = 0.18 we 
see in Fig. QJb) that all curves overlap and result 
therefore in the same Tb value, as expected. Analo- 
gously, for interacting conditions and under the same 
H/Ha — 0.18 value, the ZFC curves separate if us- 
ing the same c — 0.10 value and so stand for differ- 
ent Tb values; however, if using the same c/cq = 0.1 
value the curves overlap and Tb is the same. We 
have therefore proven the suitability of cq as a scal- 



ing factor for the dipolar interaction analogous to 
Ha for the magnetic field. 

To illustrate the inner scale of energies as the ori- 
gin of this scaling behavior, we show in Fig. [2] the 
total energies corresponding to the curves displayed 
in Fig. [TJ recorded during the ZFC processes while 
heating the system up. 
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FIG. 2. (Color online) ct = Et/2KV data correspond- 
ing to the ZFC curves displayed in Fig. [T] 

In Fig. [2j we observe that the energy curves of the 
non-interacting case differ if using the same absolute 
field (a), but overlap if using the same reduced value 
H/Ha = 0.18. Analogously, for the interacting case 
(and same H/Ha — 0.18), the curves differ if us- 
ing the same absolute concentration (c) and overlap 
if using the same normalized value c/cq. This re- 
sult demonstrates the existence of an inner scale of 
energies as a function of Ha and Co . 

We focus now on the Tb (ts in reduced units) vs. 
c problem. Following the same procedure as in Fig. 
[U we simulate ZFC processes for the different cases 
of TableUand systematically vary c, evaluating ts as 
the maxima of the curves. In Fig. [3][a), we observe 
that the ts(c) data follows different trends under 
the same absolute values H = 100 Oe and c = 0.10, 
but that share however a common origin in the non- 
interacting c = 0.0 case if using the same reduced 
field H/H A = 0.18 (Fig. EJb)). It is worth to em- 
phasize the different trends observed in the tB vs. c 
data, which reproduce a rich variety of tendencies as 
mentioned in the introduction. A complete charac- 
terization of all different types would need however 
a more detailed analysis of the specific shapes of the 
curves, however out of the scope of this work. 

In Fig. |3jc) it is shown that the apparently dif- 
ferent trends of Fig. [DJa) share however a more 
similar tendency if plotted vs. c/cq, even for the 
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FIG. 3. (Color online) t B vs. c for, (a) H = 100 Oe; (b) 
H/H A = 0.18; and vs. c/c for, (c) H = 100 Oe; (d) 
H /Ha = 0.18, for the different particle types of Table U 



same absolute field H = 100 Oe. Furthermore, 
all of them pretty much overlap if using both the 
same H/Ha = 0.18 value and plotted vs. c/c , as 
shown in Fig. E^d). These results suggest (within 
the precision of the results) the existence of a gen- 
eral ts{H I Ha] c/co) curve, with the different trends 
of Figs. 0Ha),(b) being only different fractions of it. 

It might be argued, however, that the precision 
of the ts (c) data is not enough to demonstrate the 
scaling of the SPM phenomena, since there is some 
deviation between the is(c/co) curves. To erase 
such possible doubts we have followed a different 
approach: we analyzed the maxima of the curves 
at ts, m(<s), which is far more precise to evaluate 
than ts- The results displayed in Fig. |4] demon- 
strate with high precision the role of the co param- 
eter as a scaling factor for the magnetic dipolar in- 
teraction: the mits) data, which follows different 
trends if using absolute H — 100 Oe values, as 
shown in Figs. H{a),(b); or absolute c = 0.1 val- 
ues, as in Figs. |U[a),(c); however perfectly overlaps 
in a common m(ts) curve if using the same relative 
H/Ha = 0.18 value and plotted vs. c/cq, as shown 
in Fig. Hd). 

In summary, we report the existence of a general 
scaling factor for the influence of the dipolar inter- 
action energy in single-domain entities. It is defined 
by the dimensionless parameter cq = 2K/Mg, and 
together with the anisotropy field Ha serves to de- 
fine a general behavior for such systems as governed 
by the dipolar and Zeeman energies. Specifically, we 
applied this parameter to study the controversial re- 
sults on the evolution of Tb vs. c in assemblies of 
SPM nanoparticles, finding that the apparently dis- 



FIG. 4. (Color online) (a) Magnetization at ts vs. c 
for different co-values at a fixed field H = 100 Oe; (b) 
same data as in (a), plotted vs. c/cq; (c) same data as 
in (a) but using the same reduced field h = 0.18 for the 
different Co cases; (d) same values as in (c), plotted vs. 

c/c . 



similar tendencies reported in the literature, belong 
in fact to the same general Tb(c/cq) curve. It re- 
mains as the subject of a future work the detailed 
analysis of the shape of such Tb{c/cq) curve. 
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